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ABSTRACT 

We present a formulation of boundary condition for flows 
with small disturbances. We test our methodology in an 
axisymmetric jet flow calculation,, using both the Navier- 
Stokes and Euler equations. Solutions in the far field are 
assumed to be oscillatory. If the oscillatory disturbances 
are small, the growth of the solution variables can be pre- 
dicted by linear theory. We use the eigenfunctions of the 
linear theory explicitly in our formulation of the boundary 
conditions. This guarantees correct solutions at the bound- 
ary in the limit where the predictions of linear theory are 
valid. 

Keywords: Nonreflecting boundary condition, Com- 

putational aeroacoustics, Jet flow computations 


INTRODUCTION 

Any attempt to directly compute the noise source from 
the flow field demands high accuracy of the numerical meth- 
ods, including the treatment of the boundary conditions. 
Treatment of the outflow boundary for stable and accurate 
flow simulations has attracted considerable attention [see 
for example Engquist and Majda (1977), Bayliss and Turkel 
(1982), Scott and Hankey (1985), Hagstrom and Hariharan 
(1988), Roe (1989), Giles (1990), Hariharan and Hagstrom 
(1990), Thompson (1990), Tam and Webb (1993), Atkins 


and Casper (1994)]. One usually idealizes a physical prob- 
lem to formulate the conditions at the boundary. The ef- 
fectiveness of the boundary condition is dependent on the 
degree of validity of the idealized assumptions in the actual 
flow situation. Approaches based on linear analysis, espe- 
cially variations of the characteristic methods, are widely 
used. Various investigators have derived essentially the 
same asymptotic pressure boundary conditions [See Hay- 
der and Turkel(1994) for a discussion of various works on 
this boundary condition, an evaluation of its effectiveness, 
and comparisons with other boundary conditions]. Hayder 
and Turkel (1994) observed that the asymptotic pressure 
boundary condition gave reasonable results. They how- 
ever recommended a small exit region beyond the region 
of interest. Their experiments indicated that Giles(1990) 
and characteristic boundary conditions with a larger exit 
layer also yields reasonable solutions. Because of the differ- 
ence in the asymptotic forms of wave equations in two and 
three dimensions, this boundary condition is slightly differ- 
ent in three dimensions from two dimensions [see Hayder 
and Turkel (1994)]. Hariharan and Hagstrom (1990) formu- 
lated higher order forms of this boundary condition. As we 
stated earlier, a boundary condition will give satisfactory re- 
sults if the assumptions used to derive the condition closely 
follow the actual flow situation. 

In this paper we present a new approach to boundary 
treatment based on the linear stability theory. The govern- 
ing equations of the fluid flow are nonlinear. However, if 
a mean flow is excited by a small disturbance, the linear 
theory can be used to predict its growth. Also, the eigen- 
functions given by the linear theory describe the profiles of 
the disturbances after an initial adjustment region. This 
phenomenon motivates our present effort to find a bound- 
ary condition for a flow with small disturbances. We as- 
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sume the profiles of the disturbances at the outflow can be 
approximated by the eigenfunctions predicted by the linear 
theory. The particular eigenfunctions chosen would gener- 
ally correspond to the most unstable inodes. However, any 
eigenmodes could, in principle, be used. The latter may 
be relevant for forced problems, where the excited modes 
may be determined by the forcing. The boundary condition 
that is developed here should be accurate for cases where 
the linear theory accurately describes the disturbances at 
the boundary and where these are dominated by a single, 
known mode. It may not be appropriate when the nonlin- 
earity in the flow is significant. 

In Sections 2 and 3, we give, respectively, the governing 
equations for our test problem and the derivation of our 
new boundary condition. A discussion of the basic scheme 
for our test problem is given in Section 4 and we present our 
results in Section 5. 


GOVERNING EQUATIONS 

We compute the flow field of an axisymmetric jet to test 
our new boundary conditions. We solve the Navier-Stokes 
equations as given below 


Qt + -P® + G r — ■ S 


where 


Q 


= J pu 



F = r 


( , pU 

pu r - r xx + p 

puv - T xr 

\ puH — ur xx - vr xr - &T X 


G — r 


( pv 

puv - T xr 

PV 2 “T rr +P 

V pvH — UT xr — VT rr - K,T r 



Q represents the solution variables, F and G are the fluxes 
in the x and r directions respectively, S is the source term 
that arises in the cylindrical polar coordinates, and Tij are 
the shear stresses. 


DERIVATION OF THE BOUNDARY CONDITION 

The governing equations for the fluid flows are the Navier 
Stokes equations. For high Reynold number flows, the vis- 
cous effect is small and one has to decide whether the bound- 
ary treatment should be based on the Euler equations or 
the Navier-Stokes equations. The difference between the 
two approaches is not just the type of boundary conditions 
but even the number of boundary conditions that need to 
be given. For inviscid subsonic flow, one boundary condi- 
tion needs to be specified at outflow corresponding to an 
incoming acoustic wave. For supersonic flow no boundary 
condition is required. For viscous flow the equations are no 
longer hyperbolic but rather incompletely parabolic. For 
inflow, four conditions need to be specified (in two dimen- 
sions), while for outflow, three conditions need to be spec- 
ified. In particular, the number of conditions does not de- 
pend on whether the local flow is subsonic or supersonic. 
Many codes use inviscid type boundary conditions. This 
is based on the assumption that the flow in the far field 
is essentially inviscid because of the high Reynolds number 
and the lack of physical boundary layers. Hayder and Turkel 
(1993) considered a framework for implementing the bound- 
ary conditions, where the formulation is of a characteristic 
type, but where viscous effects are also partially accounted 
for. We will use that framework for our present study. At 
subsonic outflow, we extrapolate three characteristic vari- 
ables from the interior and impose one boundary condition. 
This is done by solving the following set of equations. 

p t - pca t = Ri 

Pt + pcu t = i?2 (!) 

Pt - C 2 Pt = i?3 

v t = R 4 

where Ri is determined by which variables are specified and 
which are not. Whenever, the combination is not speci- 
fied, Ri is just those spatial derivatives that come from the 
Navier-Stokes equations. Thus, Ri contains viscous con- 
tributions even though the basic format is based on invis- 
cid characteristic theory. In implementing these differential 
equations we convert them to conservation variables p , m 
= pu,n = pv and E. Assuming an ideal gas we then have 

V ? 1 -f- 

Pt = { 7 “ l)CEt + — 2 — pt ” umt ~ VUt ) 


m t upt 



For subsonic outflow we calculate Rz> R± from the 
Navier-Stokes equations and set R\ as prescribed by the 
given boundary condition. For supersonic flows, all the Ri 
at the outflow boundary can be calculated from the Navier- 
Stokes equations or else by extrapolation of all the charac- 
teristic variables from the interior. 
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In this work, we assume the solution variable Q at the 
outflow behaves as 

Q(x,r,t) - Q(x,t) + Q'(x,r,t) 

where Q is the mean and Q r is the oscillatory part of the 
variable Q and 

Q f = e a<aj [Ci cos(a r x - ut) 4 C 2 sin(a r x — wt)] 

Thus at outflow, 


fp'l 


' Pr y 


f Pi y 


u' 



U r 

A coscvt + 

Ui 

B sinwt 

v f 


V r 

Vi 


<p f - 


< Vr > 


Vi > 



where u is the excitation frequency. Here, the vectors deter- 
mining the structure of the disturbances are eigenfunctions 
of the linear stability equations. They are used explicitly 
in our boundary conditions. This form for the disturbances 
should hold if the solution is well-approximated by linear 


theory. 



Ri = 

Vt - pcu t 



w ( - A p r sin ut 

4 B pi cos wt) 


- pcu)( -A u r sin 
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Pt + pcut 
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-f pcu>( -A u r sin 
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w (- kp r sin wt 
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Let A = . 

A. U) sin ut and B : 

= B u cosut. Then 


pr - pCUr Vi + P^i 1 ( 4 1 - f 1 

-p r 4- C 2 p r Pi - C 2 Pi J l B J l Rs ) 

or 

( A y | = }_[Vi- c 2 pi -Pi - pcui } ( i?2 1 
[ B J A l Pr- C 2 p r -Pr “ P<^r ) l ^3 J 


A = {<? Pi - Pi)(pCU r 4 p r ) 4 (p* + pCUi){Pr - (? Pr) 


boundary condition [equation (2)] in other frameworks. For 
example, the framework presented in Tam and Webb(1993) 
uses the linearized Euler equations. One can implement our 
condition by replacing the condition corresponding to the 
incoming acoustic wave by equation (2). One may expect 
to see some differences with the same boundary condition is 
implemented in different ways. 

We note that a similar construction could, in principle, be 
carried out at the inflow boundary. To do so, an eigenmode 
corresponding to a left-moving wave should be identified and 
its amplitude related to the outgoing characteristic variable, 
Rl . Finally, the incoming variables could be specified, again 
using the assumed form of the disturbance. We have not yet 
explored the feasibility of this approach. 


BASIC SCHEME 

We use a high order extension of the MacCormack Scheme 
due to Gottlieb and Turkel(1976). It has a predictor and a 
corrector stage and may be written as: 

The predictor step with forward differences is 

Qi = Q? + - FT) - (F? + 2 - FT + 1 )} + A ts? 

The corrector step with backward differences is 

<?i n+1 = \lQi + Qi 

+^{7(Fi - Fi -i) - (Fi - 1 - Fi_ 2 )} + A tSi] 

b/\x 

This scheme is second order in time and becomes fourth- 
order accurate in the spatial derivatives when alternated 
with symmetric variants. We define L\ as a one dimen- 
sional operator with a forward difference in the predictor 
and a backward difference in the corrector. Its symmetric 
variant L 2 uses a backward difference in the predictor and 
a forward difference in the corrector. For our computations, 
the sweeps are arranged as 

Q n+1 = L lx L lr Q n 
Q n+2 = L 2 r L 2 x Q n+1 

Further description of our implementations can be found 
in Hayder et al.(1993) and Mankbadi et al.(1994). 


Rl = A( -p r + pcu r ) + B(pi- pern ) (2) 

RESULTS 

Equation (2) is our new boundary condition and we use this We test our new boundary condition for an axisymmet- 

valuc of i?i in equation (1) for our numerical tests. We ric jet flow calculation. Details of such calculations can be 

would like to point out that one can implement our new found in Hayder et al. (1993) We note that the flow is 


3 



unstable, and hence provides a stiff test for any boundary 
condition. Here, the initial axial velocity is specified sis 

w(r) = 1[(1 + it'oo) - (1 - u^tanhi^r - 1))] 

and the corresponding temperature is given by the Buse- 
mann -Crocco integral of the energy equation: 

T(r) = T 0 + — - — M z ( 1 - u)(u - «<»)• 

Here u Q 0 = .25 and the jet center temperature is assumed 
to be equal to the outer flow temperature, i.e., To = T^. 
The jet Mach number is M=1.5 and Reynolds number based 
on the jet radius is 364,000. We excite the inflow profile at 
location r and time t as 

W(r,t) = W{r) + eRe(W'e iu>t ) 

where W ~ (p,u,v,p) T , W is the mean and W f is the 
eigenfunction of the linear stability equations correspond- 
ing to the mean flow profile which has the most rapid 
growth rate. For our numerical tests we used w = 1.08 and 
€ = 1 CT 6 . Eigenfunctions (EF) are obtained by solving the 
linear stability equations and also using our flow code. In 
our implementation of the new boundary condition, we use 
for the eigenfunctions the average of those obtained from 
the linear stability calculations and those computed from 
our flow code, i.e., 

ip 77 ) EF Linear Theory + EF code 

EF bc = — 5 

We show contours of vorticity magnitude for Navier Stokes 
computations Figure 1. Similar computations with Euler 
equations are shown in Figure 2. We used three compu- 
tational domains. They are 60, 50 and 40 radii long. All 
three computational domains are 5 radii wide in the trans- 
verse direction. The results at the outflow boundary of the 
shorter domains i.e., Xi = 40 and 50 compare well with the 
same quantities at the same locations in the long domains. 
The present boundary condition gave satisfactory results. 
Because of very high Reynolds number used in our compu- 
tations, there is virtually no difference between our solutions 
of the Euler and the Navier Stokes equations. 
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Figure la: Long Computational Domain Solution (xl— 60) 
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Figure lb: Intermediate Computational Domain Solution (xl=50) 
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Figure lc: Short Computational Domain Solution (xl — 40) 


Figure 1: Solutions of the Navier-Stokes equations 
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Figure 2a: Long Computational Domain Solution (xl— 00) 
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Figure 2b: Intermediate Computational Domain Solution (xl — 50) 
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Figure 2c: Short Computational Domain Solution (xl=40) 


Figure 2: Solutions of the Euler equations 
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